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Abstract. Number-resolving single photon detectors are essential for the implemen- 
tation of numerous innovative quantum information schemes. While several number- 
discriminating techniques have been previously presented, the Silicon Photo-Multiplier 
(SiPM) detector is a promising candidate due its rather simple integration in optical 
setups. On the other hand, the photon statistics obtained with the SiPM detector 
suffer from inaccuracies due to inherent distortions which depend on the geometrical 
properties of the SiPM. We have simulated the detection process in a SiPM detector 
and studied these distortions. We use results from the simulation in order to interpret 
experimental data and study the limits in which available models prevail. 
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Photon-number discriminating detectors have been the subject of great interest 
over the past few years. Number resolving capabilities are the key to several quantum 
state preparation schemes [IH3] as well as in the implementation and analysis of 
quantum computation schemes [3HS]. Photon- number resolution also enables the direct 
measurement of a state's photon-number statistics, from which non-classical properties 
as well as classical-to-quantum transitions can be studied and characterized |9l-|ll|. 

Standard single photon detection techniques cannot resolve the number of photons, 
and number resolution is generally obtained by combining several number-insensitive 
detectors in a cascade or an array-like structure. One of the first photon-number 
resolving detectors introduced was the Visible Light Photon Counter (VLPC) |12j . The 
VLPC has a high detection efficiency, but requires cryogenic cooling. Superconducting 
devices which provide number resolution are available but also require low working 
temperatures [131 - 11 7| . Room-temperature based solutions using standard single 
avalanche photodiode (APD) detectors offer limited number- resolution |18| . while other 
solutions based on spatial and temporal multiplexing are experimentally demanding [19- 
[22]. 

A promising approach is provided by the Silicon Photo-Multiplier (SiPM) |23| . 
The SiPM is composed of multiple silicon avalanche photo-diodes arranged on a single 
substrate. Each detection element acts as an independent avalanche photo-diode in 
Geiger-mode, which can absorb a photon and generate a confined electric discharge. 
Signals from all the detecting elements are combined to a single readout port, so that 
the intensity of the output signal is proportional to the number of impinging photons. 
The SiPM detector offers a very good number resolution, operates at room temperature 
and is easily integrated in optical setups. 

When using SiPMs, three main factors affect the measured photon statistics. These 
effects are the relatively low detection probability, determined by the quantum efficiency 
and geometrical configuration, internal noise caused by thermal excitations resulting 
in false detection, and optical crosstalk (CT) in which a photon created by carrier 
relaxation in one detection element is detected in a neighboring element [211125] . There 
exist several models [26H28] which take these phenomena into account. These models 
aid in the interpretation of experimental results and allow the reconstruction of the 
original photon-number statistics. However, these available models handle the CT 
effect in a rather limiting manner, either by limiting the CT probability, limiting the 
overall number of CT events, or limiting the number of CT stages (ignoring crosstalk 
events generated by crosstalks). Since these models were inaccurate in describing 
our experimental results, we have written a computational model which simulates the 
detection process in SiPM detectors. Using this model, we determine the conditions 
under which available models are accurate, and examine the limits imposed by their 
assumptions with regards to interpretation of experimental data. 

In the model, the detector is represented by a two-dimensional lattice where each cell 
represents an element of the detector. The impinging photons are distributed uniformly 
across the lattice. When a photon approaches a certain cell, the cell is triggered with 



Simulations of photon detection in SiPM number-resolving detectors 



3 



some probability, 77, the detection efficiency. If triggered, each of the cell's nearest 
neighbours can also be triggered with a probability e nn , the optical crosstalk probability 
These new triggered elements can continue to trigger their remaining nearest neighbours. 
This process continues until no new cells are triggered. The dead time of the detection 
elements in available SiPM detectors is longer than the photon propagation time between 
cells [23], so our model allows each element to only be triggered once in the process. 
The statistics were obtained on a 10 x 10 lattice with 100 elements, which corresponds 
to a typical commercially available SiPM detector (Hamamatsu Photonics, S10362-11- 
100U [29]). Nevertheless, the conclusions can be inferred to other sized sampled by 
translating the absolute number of detections to the fraction of occupancy (20 photons 
in a 100 pixel detector are equivalent to 80 photons in a 400 pixel detector). The samples 
are not entirely scalable in size since the number of cells along the borders does not scale 
with the number of elements. However, simulations performed on different sized samples 
showed that this has little effect on the overall results in the range of the experimental 
parameters. 

We begin our discussion with the optical crosstalk process. Figure [1] shows 
representative runs which portray the crosstalk evolution. The detector is initially 
triggered by N trg detections of the impinging photons. These triggers initiate a crosstalk 
process which results in the triggering of additional cells. The number of crosstalk events 
shown in these examples was chosen so that it corresponds to the average number of 
CTs produced over multiple runs. As expected, the number of triggered cells increases 
with the crosstalk probability, e nn . In figures [TBI and [Tdl we present an example of the 
effect the finite number of elements has on the number of generated CTs. Increasing the 
value of N trg , beyond some critical limit, does not result in additional crosstalk triggers. 
The majority of cells attempt to trigger neighbouring cells which have already been 
triggered, and do not contribute to the evolution of the process. 

The limited number of cells can cause the number of crosstalk events to decrease 
rather than increase with the addition of more initial triggers. This result is shown in 
figure [2j When the average number of triggered elements is small, the number of CTs 
grows linearly with N trg . As N trg is increased, the effective number of neighbouring 
cells, which can also be triggered, is reduced as many of these neighbouring cells have 
already been triggered. As a result, the average number of crosstalk events deviates from 
linearity and starts decreasing after reaching some peak value. Notice that this effect, 
which is caused by the finite size of the detector, comes into play even at relatively low 
numbers of detected photons. 

In order to identify the critical point where finite size effects begin to dominate the 
crosstalk process, we consider the average number of triggered elements generated by one 
initial trigger. In a detector of infinite size the crosstalk expansion is only suppressed by 
the probability for crosstalk. When the number of elements is limited, the progression is 
also suppressed due to overlaps with previously triggered cells. The average number of 
triggered elements should therefore be constant in the absence of finite size effects and 
should decrease when the finite size begins to impose a limitation. Figure [3] shows that 
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Figure 1: Representative simulation runs of the crosstalk evolution. The initially 
triggered pixels are marked with a cross and coloured grey, and the crosstalk triggered 



pixels are shown in red. (a) N t 



trg 



1, e r 
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Number of Initially Triggered Elements 

Figure 2: The average number of crosstalks events as a function of the initial number 
of triggers for different values of e nn . The number of crosstalk events increases 
monotonically with the number of initial triggers until some peak value is reached, 
after which the number begins to decrease due to finite size effects. For e nn = 1, the 
detector is saturated regardless of the number of original triggers. The solid line is the 
function N C t s = N e i ements - N trg . 
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Figure 3: The average number of triggered elements generated by a single detection. 
This value is constant and begins to decrease when the finite number of elements begins 
to impose a limitation. 



for large crosstalk probabilities, the finite size imposes an immediate restriction, whereas 
the process expands rather freely when the crosstalk probability is below e nn ~ 0.025. 

We define the critical number of initial triggers above which the size effects must 
be considered as the point where the average number of elements triggered by a single 
detection decreases by more than 10% compared to the unaffected value. In figure H] 
we show this critical value as a function of the crosstalk probability, e nn . The number 
of crosstalk-triggered cells increases with the value of e nn , and the value of the critical 
value is therefore lower. FigureH]also presents typical crosstalk values we experimentally 
measured for a range of bias voltages. The higher the bias voltage, the higher the gain 
and the crosstalk probability [24j. For the typical range of experimental values, the 
maximal number of detected photons for which finite size effects are not significant is in 
the range 20-100. These numbers are significantly higher than the number of photons 
observed in previously reported works [26l428| [30]. This explains why the modelling of 
the crosstalk effect without geometrical features of the detector was sufficient in these 
works. 

We turn now to evaluate the number of crosstalk stages involved in the process 
before it stops. The number of stages is an indication of the amount of cells triggered due 
to crosstalk generated by crosstalk events, and is an important factor in the modelling 
of the crosstalk process. The average number of CT stages is shown in figure [5j For 
small crosstalk probabilities only a small number of neighbouring elements is triggered 
due to crosstalk and the process ceases naturally after one stage even for large initial 
values of N trg - For large values of e nn , the number of crosstalks generated by other 
crosstalk events increases and we observe a rise in the number of stages. It is interesting 
to note that reference 26 had a CT probability equivalent to e nn w 0.025 and used a 
one-stage crosstalk model. From figure [51 we can see that for this value of e nn , the 
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Figure 4: The critical number of detected photons above which finite size effects become 
significant as a function of the crosstalk probability (left axis, blue circles), and typical 
crosstalk values, e nn , as we measured over a range of bias voltages (right axis, red 
squares). The measurements were conducted using a SiPM detector (Hamamatsu 
Photonics, S10362-11-100U). 



number of stages indeed does not exceed 1 even for photon numbers higher than what 
was detected in this work. On the other hand, we measured a CT value of e nn ~ 0.07 
and detected up to 14 photons. In this regime the number of stages exceeds 1. In fact, 
a simplistic 1-stage model would result in an error of over 20% in the evaluation of the 
number crosstalk events. The finite size effects are also apparent in this graph. Similar 
to the behaviour observed in figure [2l the number of crosstalk stages begins to decrease 
slowly after reaching some peak value. 

We now consider the second effect which governs the behaviour of SiPM detectors, 
the inherent loss mechanism. Apart from the low detection efficiency, when two or more 
photons approach the same detection element, only one avalanche can be generated 
and the number-resolution is lost. This non-linear response of the detector, affects 
the detection probabilities. The average number of detected photons can be generally 
written as |29] 



(.N detected) -^elements ^ 



/ rj Nph t ons \ 
1 — exp I — I 

\ *i elements J 



(1) 



where Nphotons is the number of impinging photons, t] is the detection efficiency and 
N elements is the number of elements. In the limit r) ■ N photons < N dements , where finite 
size effects can be neglected, the expression reduces to the linear relation 

{^detected) = V ' Nphotons- (2) 

We use this relation to determine a critical condition below which the finite number 
of the detection elements do not affect the detection probabilities. In figure [6] we 



Simulations of photon detection in SiPM number-resolving detectors 



7 




20 40 60 80 100 

Number of Initialy Triggered Elements 



Figure 5: The number of crosstalk stages as a function of the number of original triggers 
for various crosstalk probability values. The inset shows the number of stages for values 
of e nn in the experimental range 0.01-0.08 in steps of 0.01. 




Detection Efficiency, r| 

Figure 6: The number of detected photons as a function of the detection efficiency as 
obtained by our model. Solid lines represent the linear relation, (Ndetected) = V ^photons, 
for the respective number of impinging photons. The value of N p h otons appears above 
each plot. The dash lines are function plots of equation [TJ 

show the average number of detected photons as a function of the detection efficiency. 
We identify deviations of over 10% from a linear slope when the fraction of triggered 
elements, % photon3 ; exceeds 20%. Interestingly, the deviations from linearity depend on 

^elements 

the detection efficiency and this behaviour may be used in order to obtain the absolute 
detection efficiency of the detector. 

In order to demonstrate the use of our model for reconstructing the original photon 
statistics from measured data, we present in figure [7] thermal photon-number statistics 
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Figure 7: Measurements of thermal photon-number statistics taken using a SiPM 
detector (Hamamatsu Photonics, S10362-11-100U). The experimental data was fit to 
the statistics obtained using our computational model (solid lines), the 1-stage model of 
reference 26 (dashed line) and the recursive model of reference 28 (dot-dashed line), with 
crosstalk probabilities of (a) e nn = 0.078 ± 0.001, e\ s tage = 0.5 ± 0.3, e rec = 0.06 ± 0.02 
(green circles), and (b) e nn = 0.010 ± 0.001, e lstage = 0.1 ± 0.1, e rec = 0.02 ± 0.01 (blue 
squares). Experimental errors are smaller than their respective symbol size. 

measured using a SiPM detector. The measurements were conducted on a single 
polarization mode from a collinear type-II parametric down-conversion source |31| with 
a pulsed pump at a repetition rate of 250 kHz. The detector was operated using two 
different bias voltage values, which changed the crosstalk probability between the two 
measurements. The photon-number distribution which should appear as a straight 



crosstalk effect. Zero photons cannot generate crosstalk, and thus the probability of zero 
photons is not affected. However, the probability of measuring two or more photons 
grows considerably due to crosstalks. We fit the experimental data to the photon- 
number statistics obtained by our computational model, the one-stage crosstalk model 
of reference 26 and the recursive model of reference 28. The crosstalk probability defined 
in previously presented models is defined as the overall probability that crosstalk will be 
generated, rather than the probability of triggering a specific neighbour. We associate 
this crosstalk value with our defined e nn through the relation e = 1 — (1 — e n „) 4 , four 
being the number of nearest neighbours in the square lattice. When the crosstalk effect 
is weak, all models fit the data. When the crosstalk effect is strong, the 1-stage model 
and the recursive model show large deviations from the data. An analytical model which 
provides a good description for all CT values is currently under work. 

In conclusion, we have simulated the detection process in SiPM detectors and have 
shown that it is highly affected by the finite number of detection elements. These effects 




experiences a change in slope due to the 
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must be taken into account when the photon-number statistics of impinging photons is 
reconstructed. We have shown that simplistic modeling which does not account for the 
geometrical properties of the detection elements is applicable provided the number of 
detected photons is below some threshold which depends on the fraction of triggered 
elements and on the crosstalk probability. We have also shown the number of crosstalk 
stages which must be taken into account in order to properly model the crosstalk effect. 
Finally, we demonstrated the application of our computational model on experimental 
data. 
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